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Abstract 

In the present paper, we study the braneworld scenarios in the presence of quintom 
dark energy coupled by gravity. The first-order formalism for the bent brane (for both 
de Sitter and anti-de Sitter geometry), leads us to discuss the shape invariance method 
in the bent brane systems. So, by using the fluctuations of metric and quintom fields 
we obtain the Schrodinger equation. Then we factorize the corresponding Hamiltonian 
in terms of multiplication of the first-order differential operators. These first-order 
operators lead us to obtain the energy spectrum with the help of shape invariance 
method. 
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1 Introduction 



It is known that all analytically solvable potentials in quantum mechanics have the property 
of shape invariance [1]. In fact shape invariance is an integrability condition, however, one 
should emphasize that shape invariance is not the most general integrability condition as not 
all exactly solvable potentials seem to be shape invariance to [2, 3]. An interesting feature 
of supersymmetric quantum mechanics is that for a shape invariant system [4, 5] the entire 
spectrum can be determined algebraically without ever referring to underlying differential 
equations. 

In the present paper we would like to use this method for an interesting problem in cosmology. 
Here we consider the quintom model of dark energy [6] in the bent brane background. One 
of the most important problems of cosmology, is the problem of so-called dark energy. The 
type la supernova observations suggests that the universe is dominated by dark energy with 
negative pressure which provides the dynamical mechanism of the accelerating expansion 
of the universe [7, 8, 9]. The strength of this acceleration is presently matter of debate, 
mainly because it depends on the theoretical model implied when interpreting the data. 
Most of these models are based on dynamics of a scalar [10, 11, 12, 13] or multi-scalar fields 
[6]. Primary scalar field candidate for dark energy was quintessence scenario [10, 11], a 
fluid with the parameter of the equation of state lying in the range, — 1 < w < = ^-. The 
analysis of the properties of dark energy from recent observations mildly favor models with 
w crossing -1 in the near past. Meanwhile for the phantom model[12] of dark energy which 
has the opposite sign of the kinetic term compared with the quintessence in the Lagrangian, 
one always has w < —1. Neither the quintessence nor the phantom alone can fulfill the 
transition from w > — 1 to w < - 1 and vice versa. But one can show [6] that considering 
the combination of quintessence and phantom in a joint model, the transition can be fulfilled. 
This model, dubbed quintom, can produce a better fit to the data than more familiar models 
with w > —I. 

This paper is organised as follows. In section 2 we consider the quintom model of dark 
energy in the background of a five-dimensional space-time with warped geometry. Then we 
consider the fluctuations of the metric and quintom scalar fields. In section 3 we review the 
supersymmetry algebra with the central charge and shape invariance method. In section 
4 we obtain the factorized Hamiltonian for the bent brane, which leads us to investigate 
the shape invariance method with considering the central extended algebra. Finally by 
takeing advantage of shape invariance method we obtain energy spectrum in our interesting 
geometry. 

2 The quintom model in the bent brane background 

The quintom model of dark energy [6] is of new models proposed to explain the new astro- 
physical data, due to transition from w > — 1 to w < —1, i.e. transition from quintessence 
dominated universe to phantom dominated universe. Containing the normal scalar field (p 
and negative kinetic scalar field x, the action which describes the quintom model is expressed 
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as the following form 

S = j d A xdy^\{- l -R + \d a 4>d a 4> - \d aX d a X ~ V(<f>, X )). (1) 

where we have not considered the lagrangian density of matter field, and we take 4nG = 1. 
The line element of the five-dimensional space-time can be written 

dsl = g ab dx a dx b = e 2A ds\ - dy 2 (2) 

where a,b — 0,1,2,3,4, and e 2A is the warp factor. dS 2 represent the four-dimensional 
metric: 

ds\ = dt 2 - e 2VXt (dxj + dxj + dxj) (3) 

where A is four dimensional cosmological constant. We note that constant A is positive 
for de Sitter (dS^) spacetime, negative for anti - de Sitter (AdS^) spacetime and zero for 
Minkowski (M4) spacetime. 

At first we consider the interacting case with A = 0, also the functions A, 0, x are 

Av)My)>x(y)- 

From the Einstein and EulerLagrange equations we obtain, 

A" = -1(<P' 2 -X' 2 ), 

A' 2 = ^(0' 2 -X /2 )-^(0,x), 
= <f>" + AA'<f>', 

V x = -x"-4A' x ' (4) 
where a prime denotes a derivative with respect to y, and 

v dV v dV (*\ 

v * = W Vx = ~dx" (5) 

In order to obtain the first-order equation, we use [14] 

A> = -\ Wi 

4>' = \w h 

X' = ~\w x (6) 
From (4) and (6) the explicit form of the potential is 

V(<P,X) = 1(W 4> 2 -W X 2 )- V (7) 



Next we consider the general case with A ^ and we obtain 
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A' 2 -Ae- 2A = l(<P' 2 -x' 2 )-lv(<P,x) (8) 

The cosmological constant leads us to define the function which corresponds to the scalar 
fields and %. It means that this function is completely coupled and generally responsible 
for the cosmological constant. Thus we gain 

A' = --W--A-/Z, 

3 3 [ ' 

0' = ^ + ^A(a + 7 )^, 

X = -\w x - l -A{ 1 + (3)Z X) (9) 

where Z — Z(<f), \) is a new and arbitrary function of the scalar fields and respond for the 
presence of the cosmological constant, a, (5 and 7 are constants. The corresponding potential 
is, 

V(<f>,x) = \m + A{a + 1 )Z^)(W^ + A{ 1 -Za)Z^)- l -{W + A 1 Z) 2 

- l(W x + A((3 + 1 )Z x )(W x + A( 1 -3(3)Z x ) (10) 
Here we assume that 

= (11) 

So, by inserting this potential in the equations of motion one can obtain the following 
constraint, 

aW H Z+ + 6W*Z H + 2Aa(a + i)Z+Z H - ^aZ^W + A 1 Z) = (12) 

PW XX Z X + (3W X Z XX + 2A(3{(3 + i)Z x Z xx + pZ x (W + A 1 Z) = (13) 
For simplicity, we consider Z(<f>,x) = W((f>,x)- This possibility leads to equations, 

hw H — W — 0, (14) 



--d'W xx + W = (15) 



3 - r - 

2 

where d = 1+ ^ a) and d' = - . If we consider the combination of Eqs.(14, 15) we 

obtain: 

dW H - d'W xx = (16) 
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To explain this new result, we take super-potential as, 

^(0, X ) = 3asinh(6(/> + cx) (17) 

where b = and c = y^?, and a is a constant. By substituting super-potential (17) into 
Eqs.(lO), (9) one can obtain following equations respectively for potential, scalar fields, and 

Mv) 

V{4>, X ) = -\ a \ l + A 7) (2 + A(2 7 + 3a + 3/?)) cosh 2 (fe</> + c X ) + 3a 2 (l + A 7 ) 2 (18) 
also the 4>(y), xiv) an d Mv) as follow, 



/ 3^ 

4>{y) = ±J — In [tan(a(l + Aj)y)} 



(19) 



3d 

X (y) = ±\ — In [tan(a(l + A 7 )y)] 



X - sin(2a(l + Ai)y) 



(20) 
(21) 



Now we are going to discuss shape invariance condition. This condition help us to investigate 
the stability of system. So, first we consider the fluctuations of the metric and scalar fields. 
The perturbed metric is, 



ds 2 = e 2A (g fJ/l/ + eh lll y)dx^dx 1 ' — dy 2 

We use the coordinate z which is defined by following expression, 

dz = e~ A ^dy 

So, one can obtain the coordinate z as a follow, 

1 



and y is 



y 



a(l + A7) 
1 



a(l + A7) 



ln [tan(a(l + Aj)y)} 



arctan(e a(1+A7)z ) 



Also, A(z) be in the form, 
where for simply, = 1 and rj — (1 + A7) 



A(z) = \ ln 



-sech(rjz) 



(22) 
(23) 

(24) 

(25) 
(26) 
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3 Shape Invariance Method 

If the ground state energy is zero, we can factorize the Hamiltonian as, 

HM = B\g)B{g) (27) 

where g is (are) the real parameter (s), which give(s) us the potential, and B(g) is a first 
order differential operator. The ground state of Hi is annihilated by B{g). The partner 
Hamiltonian of Hi will be obtain with reversing the order of B and B\ 

H 2 (g) = B(g)B\g). (28) 

The spectrum of Hi and H 2 is degenerate. The only difference is that Hi has a zero-energy 
state and H 2 in general does not, so we have, 

H 2 B = BHi. (29) 

If we had for n > 0, 

H^W = E^ n W (30) 

this implies that, 

H 2 (B* n W) = E n W(B* n <»). (31) 

SO, the relation between the eigenvalues and eigenfunctions of the two hamiltonians Hi and 
H 2 are, 

*i 2 ) oc M'li, (32) 
where the ground state wavefunction for Hi(orH 2 ) can be obtained as, 

b¥q\x) = V { o\x) = Nexp(- j X W(y)d(y)) 
B^¥q\x) = ^ { o\x) = iVexp(+ j X W(y)d(y)). (33) 

We know that, supersymmetry give relationships between Hamiltonian Hi and H 2 , where 
Hi and H 2 are partner of each other. 

«=( H ( ; I)- < 34 > 

This matrix can be obtained from anticommutator H = {Q,Q\}, where Q and Q\ are 
supercharges, given by, 
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In this algebra we have 

[H,Q] = [H,Q*]=0 

{Q,Q} = {Q\Q j } = o (36) 

A potential is said to be shape invariant, if its supersymmetry partner potential has the same 
spatial dependence as the original potential, with possibly altered parameters. Suppose that, 
these Hamiltonian are linked by the condition, 

B{ 9l )B\ gi ) = B\g 2 )B{g 2 ) + c(g 2 ), (37) 

where the real parameters g\ and g 2 are related by a mapping / : gi — > g 2 , and c(g) is a 
c-number that depends on the parameter(s) of the Hamiltonian, so we have, 

H k = B\g k )B(g k ) + c{g k ) + ■■■ + c(g 2 ), (38) 

where 

9j+i = f(9j) 

B\g k )H k+1 = H k B\g k ) (39) 
The ground state of each of these sectors satisfies a first-order equation, namely 

B{g k )^ 1 {x-g k ) = 0. 



Now, we study supersymmetry with central charge. Supersymmetric quantum mechanics 
[10] can be formulated as a one - dimensional supersymmetric quantum field theory. A 
bosonic field is, then, a real - valued function of time, and a fermionic field is a Grassman - 
valued function of time. The d — 1, TV = 1 superalgebra with a central charge is specified 
by the following relations, 

{Q,Q^} = h 

[H, Q] = [H, Q f ] = 

{Q,Q} = {Q\Q ] } = c, (40) 

where Q and C are supercharge and central charge respectively. The above algebra implies 
[Q,C\ = [Qi,C\ = 

To realize the algebra (40), we represent the supercharges as matrices, 

where A is a real c-number. This approach is first to present an implementation of this 
algebra in a two - sector model, and then to generalize this construction to an arbitrary 
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number of sectors. 

The corresponding Hamiltonian for two-sector is, 

' B^B + 2\ 2 \ 

BB* + 2A 2 I ' 



H = 



C 




C> 



(42) 



To construct a model with 4 sectors , one can concentrate 2 two - sector model. It has 
supercharges 



Q 



( -Ai 


V o 




Ai 








-A 3 
5 3 



o \ 




A 3 J 





-A x 










\ 







Ai 


















-A 3 


5a f 




V 











A 3 


/ 



(43) 



we obtain, 



H 



and 



/ BjB 1 + 2A X 2 






C 





+ 2Ai 2 









BjB 3 + 2A 3 2 

B 3 Bj + 2A 3 2 



(44) 



/2A! 2 



V 











2Ai 2 








2A 3 2 








2As 2 



(45) 



As we see the sectors 1 and 2 are degenerate,with energies bounded from below by 2A 2 , 
and sectors 3 and 4 are degenerate, with energies bounded from below by 2A 2 . The only 
exceptions are that sectors 1 and 3 each have states that saturate their respective energy 
bounds while the even sectors do not. 

This suggest an enhanced algebraic structure. In four sector case,we define the shift operator 
S by 

/ 
B 1 
D 








B 3 






o / 



(46) 



We can choose D such that shape invariance and [H, S] = satisfy to corresponding case. 
For this, we suppose a unitary transformation which is represented by an operator Q such 
that B 3 = Q^BiQ. And also, we use a unitary operator U such that U 2 = Q, and the 
conserved shift operator takes the form, 



S = 



( \ 

B x 

WBiU 

y U^B X U 2 J 



(47) 
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From conservation of S, one can obtain the shape invariance relation, 

£i#i t - WBjBiU = k, 
2A 3 2 = 2Ai 2 + k + rfhU. 



(48) 



Therefore, the Hamiltonian of the four sector model is related to S in an especially simple 
way. In particular 

H = S j S + F, 



where 



/2A X 2 



V o 




2Ai 2 + A; 

2Ai 2 + k + WkU 








#4 



\ 



(49) 



(50) 



In the first three sectors, the energies are constrained by a Bogmol'nyi bound, Hj, > (F)^, 
because each of the first sector has to be degenerate with the Bogmol'nyi - saturating ground 
state of one of the first three sectors. The constants in F represent not only the Bogmol'nyi 
bounds of the various sectors, but also the first three energy eigenvalues of the original 
Hamiltonian. 

Next section, we take advantage from above information, and also apply to bent brane with 
quintom dark energy. 



4 The Stability of System with Shape Invariance Method 



In order to discuss the shape invariance condition we have to make the corresponding 
Schrodinger equation in terms of coordinate z 



d 2 jj(z) 
dz 2 



+ V(z)ip(z) = n 2 ^{z) 



where 



V{z) = -\\+\a»{z) + \a*{z) 



Now we can factorize the corresponding schrodinger equation in the form 



d 3 ... . 
T z + 2 



■i + l A '^ 



4>(z) = (n 2 + -A)4>{z 



The corresponding potential in term of z is , 

V(z) = s 2 rf — s(s + l)r] 2 sech 2 (r]z) 

where s — |. 

By using the new variable x = \z the Schrodinger equation can be written as, 

,2 „/„ , i\„„„l2/ 



d 2 ib(x) r 9 , , ,,, .1 , r n 2 9A , . , 
+ s 2 -s(s + l)sech 2 {x) tp{x) = (— + -^)i/>(x) 



dx l 



(51) 
(52) 

(53) 

(54) 



(55) 



and we have 



H 



d 2 

— s(s + 1) sec h 2 (x) + s 2 



dx 2 



Now we are going to factorize H in terms of lowering and raising operators, respectiv* 



B = — — stanh(x) 
dx 

= stanh(x), 

dx 



and one can obtain the paired Hamiltonians 

d 2 



E x = B ] B = 
H 2 = BB ] = 



dx 2 
dx 2 



s(s + 1) sec h 2 (x) + s 2 
s(s — 1) sec h 2 (x) + s 2 , 



where 

H 2 (s) = H 1 (s-i) + c(s). 

This relation shows us there is a shape invariance condition with c(s) = 2s — 1. 
In the case of a central charge, we choose unitary operator U as follows 

C/ = exp(|), t/*=exp(-|), 



where 

From Eqs. (46) and (47) we have, 

S^S 



rff( 8 )U — > f(s - 1) 














\ 





Wb^b^ 























I 








#4 


J 



with 



B 1 j B 1 
U^B^B X U 



d 2 



dx 2 
dx 2 
dx 2 



s(s + l)sech 2 (x) + (s - l) 2 

(g-l)(g) sech 2 (x) + (s-2) 2 

(s - 2)(s - 1) sec /i 2 (x) + {s- 3) 2 . 



/ 9A 



Also, by using Eqs. (49), (50), one can obtain F as follows: 

9A 

-M +2s -i o 


V o o o 



F = 



\ 




-J^ + 4s-4 



H 4 J 



10 



Therefore, the energy spectrum of Hi is 



9A, 



, 9A 1 9 . 
^ = + 

/ 9A 9n 

£r = (- T -?/) 

4 1} = H 4 (64) 

So, here we can discuss three cases as A = 0, A < and A > which are corresponding to 
flat, AdS and dS space respectively. The energy spectrum of Hi for flat, AdS and dS spaces 
are as following respectively. In the another term for A = we have 



^0 


= 






1 2 


> 




= -v 2 


< 




= H 4 





(65) 



for A < we have 



& = (-^)>o 

E { 2 ] = (-^ - (1 + A 7 ) 2 ) < Oor > 

= H 4 >0 (66) 

and finally for A > we have 

& = (-^) < o 

^ = (-^ + ^(1 + A 7 ) 2 )<0or>0 

= (-^-. 2 )<o 

EP = H 4 >0 (67) 

We note that in the all above cases we have some tachyonic states with negative energy. 
From (66) one can see, if |^| > (1 + A7) 2 then for the AdS 4 case, all states have positive 
eigenvalues. In this case the transition from AdS 4 to M 4 and dS 4 geometry is not stable. 
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5 Conclusion 



In the present paper we have described the algebra which gives a natural framework for 
understanding the origins of shape invariance in our interesting problem. The study of 
shape invariance solutions can be done by the factorization method. Our aim was to solve 
and discuss the stability of a bent brane in the presence of quintom dark energy and in 
different geometries with a non-zero cosmo logical constant. We have done the perturbation 
to the metric and fields and achieved the corresponding Schrodinger equation which was the 
second-order equation. Then we have factorized the equation to the first-order equations 
which are raising and lowering operators and have generated the algebra. From first order 
equations we easily discussed the energy spectrum and also the stability of the system in the 
transition to different geometries. 
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